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Abstract 

The aim of this contribution is to study the particle dynamics in 
a storage ring under the influence of noise. Some simplified stochastic 
beam dynamics problems are treated by solving the corresponding 
Fokker-Planck equations numerically. 

1 Introduction 

Colliders have become an important tool in high energy physics. For ex- 
ample, HERA, the electron-proton collider at DESY, consists of a 820 GeV 
proton ring and a 30 GeV electron ring. At the intersections of these two 
rings the colliding beam experiments HI and ZEUS are located. Besides these 
colliding beam experiments there are also two internal target experiments: 
HERA-B probing the proton beam halo with a wire target and HERMES 
using the longitudinally polarized electron beam. The experiments require 
optimal performance of the collider i.e. maximum luminosity (collision rates), 
high polarization degree and controled beam halo. In order to achieve these 
requirements one needs a good understanding of all phenomena and effects 
which cause a degradation of the beam quality. The beam constitutes a com- 
plicated many particle system of about 10 13 charged, ultrarelativistic (v m c) 
particles with spin which are distributed in 180 bunches. This ensemble is 
subject to external electromagnetic fields (dipoles, quadrupoles, multipoles 
and rf fields), to space charge fields and wakefields. Furthermore, various 
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scattering mechanisms (restgas, intrabeam) must be taken into account, and 
in the case of the electrons, radiation phenomena must be included. In ad- 
dition to these effects there are also various sources of noise such as rf noise, 
power supply noise, random ground motion and quantum fluctuations due 
to the radiation. Altogether, the beam constitutes a complicated nonlinear, 
explicitly stochastic many particle system. 

The goal of accelerator physics is to describe and understand via suitable 
models the dynamical behaviour of such a system - in the ideal case in terms 
of macroscopic dynamical variables such as particle density and polarization 
density (in phase space or configuration space). Mathematically, the system 
can be described via a stochastic Liouville equation (for a recent discussion 
of this approach in accelerator physics see |TJ) or via a stochastic differential 
equation. 

In the following we will restrict our considerations to the latter case, the 
Langevin-like description of dynamical systems. For more information about 
the conceptual foundations of these statistical dynamics problems we refer 
the reader to |2|, |3|. 

Our paper is organized as follows: In section 2 we summarize some ba- 
sic facts about stochastic beam dynamics in accelerators and we remind the 
reader of some mathematical results concerning stochastic differential equa- 
tions and the Fokker-Planck equation. In section 3 we apply the Fokker- 
Planck description to certain (simplified) accelerator problems and models 
such as beam-beam interaction, rf noise and diffusion out of a stable rf- 
bucket. Section 4 consists of a summary and a list of open problems for 
future work. 



2 Stochastic dynamics 

In the classical approximation and for the problems to be studied below, 
the particle dynamics in accelerators can be written in the form of a multi- 
plicative stochastic differential equation with a Gaussian white noise vector 
process £(t) 

j t y(t) = f(y,t) + T(y,m) (i) 

where y(t) = (x(t),rj(t)) consists of the n-dimensional phase space vector 
x(t) (n=2,4 or 6 for the orbital motion and n=8 for the spin-orbit mo- 
tion ||) and an m-dimensional Ornstein- Uhlenbeck type stochastic process 
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T](t), f(y,t) is an (n+m)-dimensional known vector function and T(y,t) is 
an (n+m) x (n+m) matrix (for more details see for example the review [[J). 
In the case of protons, equation (|1|) describes the stochastic Hamiltonian 
dynamics of the coupled synchro-betatron oscillations, and in the case of 
electrons - because of radiation phenomena (radiation damping and quan- 
tum excitation) - (P describes a stochastically and dissipatively perturbed 
Hamiltonian system ||. 

With the Ito and Stratonovich calculus one has the mathematical tools to 
study these multiplicative stochastic differential equations, whose solutions 
are in general Markovian diffusion processes [|7|, || |) H0> HH- ^ n this case, 
instead of studying the stochastic differential equations (|l]) directly [12|, [13 



one can also study the corresponding Fokker-Planck equation. The Fokker- 
Planck equation is a partial differential equation for the probability density 
p(y, t) and the transition density p(y, t\y Q , t ) of the stochastic process defined 
by equation ([I]), and in the Ito interpretation of this equation it takes the 
following form 

ajj>(y.*) = -E g^lMv,t)-p(y,t)]+ 

If one integrates p(y,t) over the Ornstein-Uhlenbeck-type variables rj, one 
obtains the probability p(x,t) = J p(x,r),t)dr) of finding the system at time 
t between x and x + dx in phase space. Using Np(x,t)dx = dn(x,t) where 
N is the total number of particles and where dn(x, t) denotes the number 
of particles in the volume element dx, p(x, t) can be interpreted (up to a 
constant) as phase space density of system ([!]). 

(HD is a partial differential equation with (n+m+1) independent variables 
and a solution of this equation requires initial conditions and suitable bound- 
ary conditions (natural (±00), periodic, reflecting or absorbing). Only few 
exact analytical solutions are available (mainly for the low (l+l)-dimensional 
case). A detailed and comprehensive review of the Fokker-Planck equation 
and a list of methods for its solution is given in [H] to which the reader is 
referred. 

In the following we will concentrate on the numerical solution of this equa- 
tion by direct discretization of space and time. The structure of the Fokker- 



Planck equation suggests an operator splitting algorithm [15, 16] which will 
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be applied in the next section to study some simplified (lower dimensional) 
stochastic problems and models in colliders. 




Figure 1: Density distribution, Q x = 0.7, near fourth order resonance 



3 Models of stochastic beam dynamics in ac- 
celerators 

In this section we consider three stochastic problems in storage rings. As a 
first example we study how an electron is influenced by the strong nonlinear 
fields of a counter rotating particle bunch (weak-strong beam-beam interac- 



tion model [[L7| , [13|, [18|). In this case, the equation of motion for the horizontal 



betatron oscillations is given by the following stochastic differential equation 
x + tx + Q 2 x x + f(x, t) = V2D£(t) 
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Figure 2: Contour plot of density distribu- 
tion 



where r is the radiation damping time, Q x is the horizontal tune, f(x,t) = 

SftCbb' 1 — —'$p(t) is the beam-beam force with the beam-beam parameter 
(bb, and is the Gaussian white noise process of strength D. 8 p (t) denotes 
a strongly localized periodic function. 

The numerical solution of the corresponding Fokker-Planck equation near 
a fourth order resonance is depicted in Fig. [T], and Fig. |2| shows a contour 
plot of the density, which shows the nonlinear resonance characteristics of 
the underlying Hamiltonian dynamics |L9 |. As initial conditions we have 
used a Gaussian distribution localized at the origin of the two-dimensional 
phase space (x, x = v). We also want to mention that this model has been 
used to compare various numerical tools to study stochastic systems such as 
cell-mapping methods, Monte-Carlo methods and finite differences pOf| . 

An important problem in proton storage rings is to study the influence 



of rf noise on the particle stability (see for example [21. |22|. 53, pi). Here, 



we investigate the diffusion out of a stable rf bucket under the influence of 
random energy losses (for example due to scattering) and under Gaussian 
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white noise. The dynamics is governed by 

4> + Vy sin(» + V 4 sin(4</>) + A p (t) = V2D£(t) 

where Vi, V4 denote the voltages of a two-rf system, and where A p (t) denotes 
the average energy loss due to scattering. 




Figure 3: Deterministic longitudinal phase space and bucket structure 

The unperturbed bucket structure is shown in Fig. ^ and Fig. |4] and Fig. [| 
show how the distribution diffuses outward in phase space thus forming a 
coasting beam. The initial distribution was Gaussian and was localized near 
the origin. 

The final example we want to show treats the longitudinal particle motion 
under the influence of coloured rf (amplitude) noise with the equation of 
motion for the phase 

4> + n 2 a (l + ri(t))mn(<t>) = Q 

where rj(t) denotes an Ornstein-Uhlenbeck process 

fj = -ar} + V2D£(t). 

D and a are parameters which define the correlation time of the process. The 
numerical solution of the corresponding Fokker-Planck equation integrated 
over the Ornstein-Uhlenbeck variable (rj) is depicted in Fig |^ for an initial 
distribution localized near the unperturbed separatrix of the system. 

More examples which illustrate the usefulness of the Fokker-Planck de- 
scription of stochastic problems in accelerators can be found in |TJ| 
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Figure 4: Density after 80000 turns, strong kick 10 5 every N=100 turns, 
and strong noise D = le — 13. 

4 Summary and discussion 

In this contribution we have shown that stochastic beam dynamics is an 
important issue in accelerator physics. Usually - via various diffusion mech- 
anisms - noise can lead to a degradation of the beam quality in a collider 
(emittance growth, reduced lifetime etc). However, the application of clev- 
erly chosen noise in the transverse plane or longitudinal phase plane can also 
help to shape and control the beam and its halo [JS^, |25|| . Furthermore, we 
have shown that the Fokker-Planck equation is a suitable and helpful math- 
ematical tool to treat these stochastic systems. Since only few exact results 
are available for this partial differential equation (especially in higher dimen- 
sions) one needs a powerful, reliable, accurate and fast numerical solver. Such 
a solver which is based on the operator splitting method has been developed, 
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Figure 5: Density after 800000 turns, weak kick 10 6 every N=100 turns, 
and weak noise D = le — 14. 



and it has been used to study stochastic beam dynamics problems in ac- 
celerators. However, longtime calculations and higher dimensional problems 
such as stochastic spin-orbit motion in realistic colliders with at least (8+1) 
independent variables are certainly beyond the capacity of this and other 
codes even if parallel algorithms and high performance computers are used. 
So -in addition- to these numerical studies of the Fokker-Planck equation 
one also needs perturbative methods such as averaging with all its mathe- 
matical subtleties. Furthermore, complementary studies of the dynamics via 
direct analysis of the underlying stochastic differential equations or discrete 
stochastic maps are important. Another alternative could be the use of ana- 
logue computers (see the recent interesting review [p6fl ). In addition to these 
topics noise induced transitions (stabilization by noise) [^TJ and stochastic 
resonance ]27| might play an important role in accelerator physics. First 
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Figure 6: longitudinal phase 
space distribution in case of 
coloured amplitude noise inte- 
grated over Ornstein-Uhlenbeck 
variable 



steps to investigate these phenomena in storage rings have been undertaken 

m 
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